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2 ■ ABSTRACT 

-*— > 

BVRIJHK photometry, Spitzer-GIAMPSY, photometry and HK band 
^ ' spectroscopy were used to study the steUar content of IRAS 19343+2026, a 

cn ' 

\j^ \ (proto) star/cluster candidate, located close to the Galactic plane. The data 

m ' 



% 



suggest that IRAS 19343+2026 is a rich cluster associated with a massive 



Xf^ i protostar of 7.6 Mq with an age of ^ 10^ yr. Three point sources in the 

^— V ' vicinity of the far-infrared (FIR) peak are also found to be early B type 

stars. The remaining (predominantly low mass) members of the cluster are 
best represented by a 1 - 3 Myr pre-main-sequence (PMS) population. HK 
band spectra of two bright and five faint point sources in the cluster confirm 
that the results obtained from the photometry are good representations of 
their young stellar object (YSO) nature. Thus, IRAS 19343+2026 is a young 
cluster with at least four early B-type stars classified as young (10'' - 10^ yr), 
that are surrounded by a somewhat older (1-3 Myr) population of low mass 
YSOs. Together, these results argue for a scenario in which low mass stars 
form prior to massive stars in a cluster forming environment. We compute 
the Initial Mass Function (IMF) for this cluster using the X-band luminosity 
function; the slope of the IMF is shallower than predicted by the Salpeter's 
mass function. The cluster mass, Mtotai, is estimated to be in the range '^ 307 
M0 (from the data completeness limit) - 585 M© (extrapolated down to the 
brown dwarf limit, assuming a certain IMF). 
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1 INTRODUCTION 

One of the important issues in understanding massive star formation i s the sequence of for- 



matio n of low and high mass stars in a cluster forming environment (jZinnecker fc Beuther 



20081 ). This issue has relevance when explaining the nature of massive star formation in 
terms of either: (a) low to intermediate mass stars accumulating matter while evolving to- 
ward the main sequence, or (b) competitive accretion of massive cor es in a star forming cloud 



flMcKee fc Tanll2002l : McKee fc Tanll2003 ; 



Zinnecker fc Yorke 



20071 ). Observations of embed- 



ded clusters generally display mass segregated configurations, with a few young high mass 
stars located at the cluster centers (JMcNamara fc Sekiguchilll986l . and references therein). 
Although an age estimation for the lower mass pre-main-sequence (PMS) population is pos- 
sible using evolutionary tracks and colour-magnitude (CM) diagrams, the same is not true 
for massive stars, since they are usually only a few in number. Therefore, it is necessary 
to investigate clusters around high mass protostellar object candidates, where the youth of 
the massive stars is made clear by vario us other signposts, such a s out flows, compact HII 
regions, excess FIR emission, and so on. iKumar. Keto fc ClerkinI ( l2006l ) searched for clus- 
tering around a sample of high mass protostellar objects and found 54 embedded clusters 
among 215 candidates. The detection of clustering in about 25% of the sources (all targets 
located away from the Galactic mid-plane) suggests that at least one generation of low mass 
stars has already formed around these massive protostellar candidates. However, the Two 
Micron All Sky Survey (2MASS) data used by these authors did not effectively uncover the 
lower mass population; the 2MASS data were not deep enough to ch aracterise the remaining 
cluster members. In contrast, a follo w-up of the same data sets by iKumar fc Gravel (120071 ) 
using the Spitzer-GLIMPSE survey (jChurchwell et al.ll2009l ) did not find any clustering for 
the sources in the Galactic mid-plane. This mid-infrared survey should have been more 
successful at finding clusters because of the better dust penetration. 

In this paper, using high quality optical and infrared observations of a high mass proto- 
stellar object in the Galactic mid-plane, we attempt to examine such discrepancies, measure 
the stellar content, and understand the sequence of low and high mass star formation in 
the associated cluster. The high quality NIR photometric and spectroscopic data are com- 
bined with Spitzer IRAC and MIPS observations to study the target. Study of the massive 
star content is particularly aided by radiative transfer modelling of the SED. The target. 
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IRAS 19343+2026, is a high mass protostellar object candidate located at a kinematic dis- 
tanc e (measured using N H.s hne velocities) of 4.2 kpc, with a FIR luminosity of 2.7x10^ 



Lf7^ flMolinar 



Mercer et al. 



et al. 



19961 ). This cluster is detected in the GLIMPSE survey (cluster 24 of 



20051 ). Section 2 describes the detailed observations. The photometric, spec- 



troscopic and SED modelling results are also presented in this section. The implications of 
the results are discussed in Section 3. We then summarize our conclusions in Section 4. 



2 OBSERVATIONS AND DATA REDUCTIONS 



2.1 J,H,K imaging and photometry 

NIR photometric imaging observations were made at 

Telescope (UKIRT) with the facility imager UFTI (JRoche et al. 



:he 3.8 m United Kingdom Infrared 



2002). UFTI houses a 



HAWAII-1 1024 X 1024 pixel array. The UFTI plate scale of 0.091" gives an available field 
of view (FOV) of ~ 90". Photometric observations through J (A = 1.25 /im, AA = 0.16 /im), 
H {X = 1.64 /im, AA = 0.29 /xm) and K {X = 2.20 /xm, AA = 0.34 /im) broad-band filters 
were obtained for the IRAS source during the night of 26*^ June 2002. An integration time 
of 60 sec was used in each of the J, H and K band filters; averaging jittered exposures 
yielded a total exposure time of 540 seconds in each band. The mean seeing measured was 
0.6" in the K-hand images. A nine point (3x3) jittered observing sequence was executed to 
obtain data that provided final mosaics with a total FOV of ~115"xll5". We note that the 
signal-to-noise ratio at the edges of these mosaics is lower than within the central 90" area. 
Standard data reduction techniques involving dark subtraction and median-sky-fiat-fielding 
of the jittered object frames were applied. The K-hand image of IRAS 19343-1-2026 is shown 
in Fig. 1, with an overlay of Spitzer MIPS 24 /im contours. 

Subsequent to our obser vations, this region wa s recently covered by the UKIRT Infrared 
Deep Sky Survey (UKIDSS: iLawrence et al.ll2007l ) Gala ctic Plane Survey (GPS). The GPS 



is an ambitious survey of the Northern Galactic plane (JLucas et al. 



20081 ). The aim of the 



survey is to map 1800 square degrees of the plane (|b| < 5 deg) in J, H, and K to a depth 
of J ~ 20.0, H ~ 19.1, K ~ 19.0 at sub-arc second resolution. UKIDSS employs the Wide 
Field Camera (WFCAM; 



Casali et al 



20071 ) at UKIRT. WFCAM contains four Rockwell 



Hawaii-II (HgCdTe 2048x2048 pixel) arrays spaced by 94% in the focal plane. With a pixel 
scale of 0.4" the field of view of each array is 13.7'. All UKIDSS survey data are reduced 
by the Cambridge Astronomical Survey Unit (CASU) and are distributed via the WFCAM 
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Science Archive (WSA|1|). The GPS data sensitivity is comparable to our observations {K 
=18 mag), albeit with a relatively poor seeing (1") and a coarse pixel scale of 0^2^^ F ig. 2 
shows the SO*'* nearest-neighbour (NN) density map (ISchmeja. Kumar fc Ferreirall2008L and 
references therein) of the UKIDSS i^-band source counts. From this figure the cluster center 
was found to be at a2ooo = 19^^36^31*^, ^2000 = +20°33 03 . For the purpose of evaluating the 
foreground/background contamination, we have chosen a "control region" ~ 2' North of the 
cluster center at aaooo = 19'^36'^32^ ^2000 = +20°34'48" (see Fig. 2). 

We utilised tasks available in the Image Reduction and Analysis Facility (iraf) package 
for our photometric analysis. DAOFIND was used to identify sources in each image. A 
point spread function (PSF) model was computed by choosing stars of different brightness 
that were well spaced out in our images. Photometry was performed using the DAOPHOT 
package. Aperture corrections were determined by performing multi-aperture photometry 
on the PSF stars. The instrumental magnitudes were calibrated to t he absolute scale usin g 



observations of UKIRT faint standard stars; FS 29, FS 35 & FS 140 (iHawarden et al 



200 ih 



These standards were observed over a range in airmass (1.05 - 1.79) that was comparable 
to the target observations. The res ulting photometric data a re in the natural system of 



the Mauna Kea Consortium Filters ( ISimons &: Tokunaga 



these data, we converted magnitudes to the 



Bessell fc Brett 



20021). Fo r the purpose of plotting 



(119881 ) (hereafter BE) system. 



since the main-sequence references are in the BE system. To do this, we first converted the 
Mauna Kea syst e m to the CIT system and then to the BE system using equations given by 
Hawarden et al.l (120011 ). 

Representative sub-images consisting of stars and nebulosity were chosen to determine 
the completeness limits. Limits were established by manually adding and then detecting 
artificial stars of differing magnitudes. By determining the fraction of stars recovered in each 
magnitude bin, we have deduced 90% completeness limits of 19.3, 19.0 and 17.5 magnitudes 
in the J, H and K bands, respectively. Our observations are absolutely complete (100%) 
to the levels of 17.3, 17.2 and 16.2 magnitudes in J, H and K, respectively. Photometric 
analysis was carried out using data with photometric errors of less than 0.1 mag. Absolute 
position calibration was achieved using the coordinates of a number of stars from the 2MASS 
catalog. The astrometric accuracy of the data presented in this work is better than 0.5". 



^ http://surveys.roe.ac.uk/wsa/index.htinl 
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The sources are saturated aX K < 12. For such bright sources, 2MASS Point Source Catalog 
data were used. 

We also extracted source magnitudes from the UKIDSS images using the methods de- 
scribed above, although our UFTI images (rather than 2MASS data) were used to calibrate 
the UKIDSS data. The UKIDSS observations were used to create a density map of the region 
(see Figure 2), and to obtain photometry over a larger area than was observed with UFTI. 

2.2 HK band spectroscopy 

NIR HK band spectra of several stars i n the target field were obtai ned using the UKIRT 
1-5 /im Imager Spectrometer (UIST) (jRamsay Howat et al.l l2004l ) on the night of 21** 



July 2005. UIST has a 1024 x 1024 InSb array and a 0.12 arcsec per pixel plate scale. 
The observations were made using a 120 arcsec long slit with a width of 4 pixels. The HK 
grism was used, which allowed complete wavelength coverage from 1.395 /xm to 2.506 /im 
with a spectral resolution of R ~ 500. The exposure time per frame was set to 60 sec; 8 
exposures were obtained per slit position, resulting in a total integration time of 8 minutes. 
The 120" long slit was aligned carefully to cover multiple sources in a single position. This 
configuration and integration time provided a 5cr point source sensitivity of ~ 16 mag in the 
HK bands. The target was nodded up and down the slit in an ABBA fashion; B exposures 
were subtracted from adjacent A exposures to give sky-subtracted 2-dimensional spectral 
images. The spectral images were averaged and the positive and negative beams of stellar 
spectra were then optimally extracted. These resulting "group" spectra were corrected for 
telluric absorption and flux calibrated by division with a similarly-observed standard star 
spectrum. The G2V source HIP 102 189 was used for this purpose. (Note that the standard 
star spectrum was first divided by a black body template of appropriate temperature to 
preserve the intrinsic slopes of the calibrated target spectra.) 

The initial reduction steps (i.e. dark and flat c orrection) were perfo rmed using the 



UKIRT/Starlink data reduction pipeline ORAC-DR fICavanagh et al.ll2003l ): the extraction 



and analysis of the spectra were done with IRAF and IDL. 

2.3 Optical Imaging and Photometry 

Bessell 5V7?/ images of the region associated with the IRAS 19343+2026 region were ob- 
tained with the 2-m Himalayan Chandra Telescope (HCT) of the Indian Astronomical Ob- 
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servatory (lAO), Hanle, India, on 16*^^ March 2005. The Hanle Faint Object Spectrograph 
Camera (HFOSC), which is equipped with a SITe 2K x 4K pixeP CCD, was used. With a 
pixel scale of 0.296", the FOV of HFOSC, where only the central 2K x 2K region is used 
for imaging, is ~ 10 x 10 arcmin^. Further details on the telescope and the instrument 
can be found at http://www.iiap.res.in/iao/hfosc.html. Observations were carried out under 
good photometric sky conditions. The typical seeing [full width at half-maximum (FWHM)] 
during the period of observations was ~ 1.8". Bias and flat frames were obtained at the 
beginning and at the end of th e observing night. Photometric standard stars around the SA 



111-775 region (JLandolt 



19921 ) were observed to obtain the transformation coefficients. 



The data reduction was again carried out using IRAF tasks. Object frames were flat- 
fielded using median-combined normalized fiat frames. Identification and photometry of 
point sources were performed using the DAOFIND and DAOPHOT tasks, respectively. 
Photometry w as obtained using the PSF algorithm ALLSTAR in the DAOPHOT package 



Stetson 



19871 ). The residuals to the photometric solution were <0.05 mag. 



2.4 Spectral energy distribution modelling 

The Spitzer-GLIMPSE survey IRAC and MIPS images and IRAC photometry of this tar- 
get were analysed using the IRSA-IPAC image cutouts and GATOR facilities. MIPS 24 
/xm photometry was extracted using the APEX single frame pipeline on the Post-Basic- 
Calibrated Data. The photometric data in the S^^zteer bands (3.6 - 24.0 /im) were combined 
with 2MASS NIR photometry and IRAS data to construct SEDs for the four brightest 
sources in th e region - the massive young star candidates. The online SED fitting tool 



developed by 



Robitaille et al. 



(120071 ) was used to fit the resulting SEDs in the 1-24 /xm 
bands. The SED fitting tool is based on matching the observed SEDs with a large grid of 
precomputed radiative transfer models. The models assume an accretion scenario with a 
star, disk, envelope and bipolar cavity, all under radiative equilibrium. While the mass and 
age of the star are uniformly sampled within the grid limits, the radius and temperature 
are interpolated using evolutionary models. The optical and NIR data points constrain the 
stellar parameters, the mid-infrared fluxes cons train the disk paramet ers and the far-infrared 



points are sensitive to envelope emission (see 



Robitaille et al. 



2006 



for more details). The 



SED fitting tool has been successfully tested on low mass stars and shown to produce reliable 
estimates of physical parameters by comparing against values obtained by other independent 
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measurements. 



Grave fc Kumar! (120091 ) discuss the implementation of this tool on a large 
sample of massive protostellar objects. For fitting purposes a 10% error was assumed on all 
fluxes. The weighted mean (weights being the inverse of x^) of each physical parameter was 
computed for all models that satisfied the criteria y^ — x^best < 3, where x^ is the statistical 
goodness of fit parameter measured per data point. 

3 RESULTS 

A rich cluster of stars can be seen in our K-hand image (see Fig. 1), the brightest stars in the 
center of the field coinciding with the peak of the MIPS 24 /xm contours. The two straight 
lines in Fig. 1 mark the slit positions used to obtain the spectra of the stars numbered 1 
to 7. Stars 1, 3 and 8 were modelled using 1-24 yum SEDs. Star 9 was not detected in the 
MIPS 24 yum band. The brightest stars in the field are expected to be massive stars, while 
the fainter population represents the low mass members. In Fig. 3 (see the online electronic 
version of this article for a colour plot) a colour composite of the Spitzer infrared images is 
shown. The three-colour composite image was made using the Spitzer IRAC 3.6 /im, IRAC 

8.0 yum and MIPS 24 yum images coded as blue, green and red, respectively. Notice the 
bipolar shape of the nebula and the FIR excess emission appearing as red. The stars 8 & 9 
are associated with excess emission at 24 yum, as can be seen from the red colour surrounding 
these two stars. Also, the rich cluster of stars seen in Fig. 1 is not visible in this composite 
image. We will comment further on this issue in Sect. 4. In the following, we first discuss the 
photometric analysis of the whole sample, then the IMF of the cluster region, followed by 
spectra of the seven representative sources (stars 1-7), and finally the SED fitting analysis 
of the brighter massive star candidates. 

3.1 Photometric Analysis 

The JHK photometry of the point sources in the IRAS 19343+2026 region was used to 
construct colour-colour (CC) and CM diagrams. A combination of CC and CM diagrams 
made using various optical and infrared bands were then used to evaluate the membership of 
the cluster against a control field, and to evaluate the cluster properties. The total number 
of point sources (with magnitude errors < 0.1 mag), detected in the region shown in Fig. 1 
was 333, 688 and 875 in the J, H and K bands, respectively. 307 stars are common to all 
three bands; 607 stars appear in the H and K bands. The UFTI JHK photometric data. 
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Figure 1. A _ft'-band image of the IRAS 19343+2026 region shown with a logarithmic greyscale stretch. Contours represent the 
MIPS 24 /im emission. The solid lines mark the positions of the sUts used to obtain HK band spectra of the stars numbered 
1—7; note that spectra were not obtained for the bright NIR sources labelled 8 and 9. The abscissa (RA) and the ordinate 
(DEC) are in J2000.0 epoch. 

along with the positions of the stars, are given in Table A.l; the complete table is available 
in electronic form as part of the online material. 

Photometry of the images obtained from the UKIDSS-GPS survey was performed over a 
larger area (~ T x T). The analysis of both the cluster and control regions was carried out 
in a similar way. In order to evaluate the mass function (see Sect. 3.2) of the cluster region, 
we have chosen a box size of 1.2' centered on the cluster and control regions, as mentioned 
earlier (see Sect. 2.1). 
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Figure 2. 50* nearest-neighbour density map of UKIDSS ft'-band source counts around the IRAS 19343+2026 region. The 
contours are plotted over the UKIDSS image; the contour level starts at 111 stars per arcmin'^ and increases in steps of 28 
stars per arcmin^. The solid and dashed boxes represent the chosen cluster and control field regions. 



Fig. 4 shows the J — H/H — K CC diagram for the cluster field. The solid {bottom-left) 
and broken heavy curves represent the main-sequence dwarf and giant stars, respectively . 
The dotted line indicates the locus of T Tauri stars (iMeyer. Calvet fc Hillenbrandl 119971 ). 
while the dashed-dotted line represents the HAeBe locus ( Lada &: Adamslll992[ ). The three 
dashed parallel lines are the reddening vectors. We adopt a slope consistent with E{J — 
H)/E{H — K) = 1.9 (BB system), wh ich is appropriate to an interstellar reddening law of 
R = 3.12 (IWhittet fc van Bredalll98Cll ). The CC diagram was used to identify the reddened 
population that falls in the region occupied by T Tauri and HAeBe stars. Stars that lie 
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Figure 3. Three-colour composite made from Spitzer images (for a colour representation see the online electronic version of 
this paper). The MIPS 24 ^m, IRAC Ch4 (8.0 fim) and IRAC Chi (3.6 fim) images are coded red, green and blue respectively. 
The abscissa and the ordinate are in J2000.0 epoch. 



outside the region of reddened main-sequence objects (i.e. redward of the reddening hne 
drawn from the base of the main-sequence dwarf branch, that is, redward of the middle 
of the three vectors) are YSOs with intrinsic colour excesses. We shall refer to these as 
"probable cluster members" . By de-reddening the stars on the CC diagram that fall within 
the first two reddening vectors (encompassing the main-sequence and giant stars) to the 
dwarf locus, a visual extinction (Ay) for each star was calculated. The individual extinction 
values range from to 20 mag, resulting in an average extinction of Ay = 7.6 mag. 

Fig. 5 shows the UKIDSS GPS K/H — K CM diagram for the cluster and the control 
region. From the control region plot {middle panel) it is possible to clearly identify the dwarf 
and giant branches of the stellar population, i.e. sources with an H — K colour of less than 
1.0. In contrast, the cluster region {left panel) displays the dwarf branch (the group of stars 
running vertically aX H — K ^ 0.4) and a significant number of stars with H — K colours 
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Figure 4. CC diagram of the region around IRAS 19343+2026 constructed from the UFTI data. The solid (bott om-left) and 
broken heavy curves represent the main-sequenc e dwarf and giant stars, respectively , taken from iBessell fc Bretj lll988 V The 
dotted line in dicates the locus of T Tauri stars JMever. Calvet &: Hillenbrandlll997l ): the dashed-dotted curve represents the 
HAeBe locus JLada &: Adam jl 19921) . The dashed parallel lines are reddening vectors on which crosses separated by Ay = 5 mag 
have been marked. The triangles mark the brightest NIR sources in the region, while the numbers refer to the sources labelled 
in Fig. 1. The star symbol (star 1) marks the 2 /im source that coincides with the FIR peak. 



around 1. These stars are essentially cluster members which are also identified on the CC 
diagram as probable cluster members. The control field also has very few stars with H — K 
< 1 as compared to the cluster region. Therefore, we find that the contamination from 
foreground/background stars in our cluster region is minimal compared to the statistics of 
the probable cluster members. However, H — K colour alone is no good for distinguishing 
between a reddened background star and a young star. Therefore, we have used a statistical 
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Figure 5. K/H — K CM diagrams of the IRAS 19343+2026 cluster region (left), the control region (middle) and the sta- 
tistically cleaned cluster region (right) from the UKIDSS GPS data. The solid and dotted curves show the PMS tracks from 
iPalla & Stahler (1999) and Baraffe et al. (1998), respectively, fo r an assumed age of 3 My r for low mass range (M < 3 Mq), 
while the dashed curve denotes the PMS track for 3 Myr age l|Leieune fc Schaerenl200ll ) for intermediate mass range (4 < 
M/Mq < 8.5). The parallel slanting lines identify the reddening vectors. 



approach to remove contaminating field stars from tlie K/H — K CM diagram of the cluster 
region (we will discuss this further in detail in the following Sect. 3.2). The cluster member 
stars are clearly seen in the statistically cleaned CM diagram [right panel). 

Our aim here is to estimate the approximate age of the lower mass population by com- 
paring our data with theoretical evolutionary tracks. The low mass evolutionary tracks have 
been tested on various observational data of young stellar clusters and are thought to be 
good representations. The evolutionary models for high mass stars are more complicated and 
varied, however, and are not well tested against observational data sets. Therefore, we first 
use the pre-main-sequence (PMS) isochrone fitting for the lower mass population. Further, 
it should be noted that a CM diagram that uses only one colour may not be very effective 
in describing the less understood massive young stars. 

Since the stellar photosphere's are better represented by shorter wavelength data such 
as the optical, J and H bands, we use these data when comparing source photometry with 
the model isochrones. In the region of Fig. 1, all optically detected stars in the BVRI bands 
were plotted on a V^ vs V^ — / CM diagram. All but seven of the optically detected stars are 
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foreground stars. Since seven sources is not statistically significant for the young cluster, we 
instead plot the J vs J — H CM diagram in Fig. 6. In this plot the slanting arrow denotes 
the reddening vector associated with 15 mag of extinction. YSOs identified in the JHK CC 
diagram in Fig. 4 are shown a s open circles; triangle s mark the bright stars. 

PMS evolutionary tracks (jPalla fc Stahlerlll999l ) for ages 1, 3 and 5 Myrs are shown in 
Fig. 6 with dotted, short-dashed, and l ong-dashed lines, respectively. We have assumed a 
distance of 4.2 kpc (jMolinari et al.lll996l ) and have reddened the isochrones with the average 
extinction of A^ = 7.6 mag estimated from the JHK CC diagram analysis in Sect. 3.1. 
Given the poor statistics for this distant cluster, it is not easy to derive the age by isochrone 
fitting. Even for regions with good statistics, it is quite difficult to constrain the age because 
most low mass stars spend the bulk of their PMS time on the Hayashi track, which is the 
verticalpart of each isochrone. However, a small fraction of observed young stars will coincide 
with the brief evolutionary phase associated with the Henyey track, which is the horizontal 
transition between the Hayashi track and the Zero- Age-Main Sequence (ZAMS) (see e.g. 



Ascenso et al. 



20071 ). 



In Fig. 6 we note that there is a large population of low mass/low J magnitude candidate 
young stars (open circles) that lie to the right of the 1 Myr isochrone (dotted line). Some 
of these sources may be Class I protostars associated with the cluster. However, most will 
be Class II sources (T Tauri stars) with an age of about 1 Myr; many of these sources lie 
to the right of the 1 Myr isochrone because of extinction which, to these embedded sources, 
will be higher than the Ay ~ 8 mags used to deredden the PMS isochrones. However, there 
is also a population of candidate low mass young stars that lie to the left of the 1 Myr 
isochrone. These sources represent a more evolved group of young stars that are probably 
3 Myr or older. Indeed, there are five probable cluster members that align horizontally and 
may therefore represent the Henyey part of the 3 Myr isochrone (short dash). Hence, we 
estimate that there is a low mass population of stars associated with IRAS 19343+2026 that 
is best represented by an age of 1 Myr or more. The mass range plotted in the figure is from 
0.6 M0 to 3.0 M0 for the 3 Myr isochrone. The 5 Myr isochrone (long dash; mass range 
0.6 - 2.5 Mq) may well be within the age spread of the cluster. However, it extends beyond 
the limits of the observed data points. 

One should also note that the J vs J — H diagram contains only 333 sources that are 
common to the J and H bands. This is a considerable under-representation of the full sample, 
given that 875 sources are detected in the K band and that many of the fainter, redder 
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Figure 6. J vs J — H CM diagram from the UFTI data. Open circles are probable cluster members as estimated from the 
J — H/H — K CC diagram in Fig. |4] The triangles mark the brightest NIR sources in the region, while the numbers refer to 
the sources labelled in Fig. 1. The star symbol (star 1) marks the 2 /im source that coincides with the FIR peak. The PM S 1 
Myr, 3 Myr (range = 0.6 Mq- 3.0 M0; marked) and 5 Myr (range = 0.6 Mq- 2.5 Mq) isochrones l lPalla fc Stahledll999l) are 
drawn using dotted, short-dashed and long-dashed lines. The slanting solid line shows the reddening vector. 

sources will be detected only at longer wavelengths. Nevertheless, the sample is statistically 
significant when comparing with evolutionary tracks and obtaining an age estimate for the 
low mass sources in this region. 



3.2 The Initial Mass Function 

As a next step to better understanding this region, we evaluate the mass function of the 
cluster. In the NIR, young stars are best detected in the ii'-band. Consequently, we use 
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the K-bandi. Luminosity function (KLF) to evaluate this statistically important parameter. 
Ideally, we want to build the KLF using UFTI data which has a higher spatial resolution 
compared to the UKIDSS data. But the estimate of the mass function should be compared 
with the mass function of the control field, which is not possible using UFTI data, since we 
do not have control field images obtained with UFTI. The number of sources detected in 
the UFTI and UKIDSS images is very similar (to within 5%) for the same FOV. Therefore, 
we will evaluate the initial mass function (IMF) using the UKIDSS data which is uniformly 
calibrated for both the cluster and control regions. 

By using the control region comparison, the differences arising due to contamination 
from the foreground and background members are accounted for. Although we identified 
probable cluster members (T Tauri and related sources) in Sec. 3.1, this identification was 
based on using only those stars that appear in the J, H and K bands. However, there are 
roughly twice as many stars detected in the i^-band, whose cluster membership can not be 
easily verified. Therefore, we have used statistical criteria to estimate the number of probable 
member stars in the cluster region. To remove contamination of field stars from the cluster 
region sample, we have statistically subtracted the contribution of field stars from the CM 
diagram of the cluster region using the following procedure. For a randomly selected star 
in the K, {H — K) CM diagram of the control region (see Fig. 5), the nearest star in the 
cluster's K, {H — K) CM diagram within K ± l.Q and {H — K) ± 0.5 of the field star was 
removed. The statistically cleaned K, {H — K) CM diagram of the cluster region is shown 
in Fig. 5 (right) which clearly shows the presence of PMS stars in the cluster. With th e 



help of this statistically cleaned CM diagram, we compute the IMF (e.g.. lOjha et al.ll2009[ ). 
The K-hand luminosities are corrected by the mean extinction determined by dereddening 
the stars on the CC diagram (see Sect. 3.1) which is Ay ~ 8.0 mag. The de-reddened 
magnitudes were then used to construct a reddening-corrected KLF. Stellar masses for the 



assumed 3 Myr old c . 



( JLejeune fc Schaerer 



uster mem bers were obtained usi ng the evolutionary models of Geneva 



20011) and 



Palla fc Stahlerl (|l999|). The MF of the statistically cleaned 



cluster region sources was obtained by counting the number of stars in various mass bins 
and is shown in Fig. 7. The vertical dashed and dot-dashed lines represent 90% and 100% 
completeness limits obtained for the UKIDSS GPS data, respectively. 

The MF of the IRAS 19343+2026 region appears to rise monotonically up to ~ 1.0 Mq. 
The MF has a slope {dLog{N) / Log{M)) of -1.12±0.16 for the probable cluster members 



sample over the mass range 1.0 < M/Mq < 6.3. The classical value derived by ISalpeter 
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fll955h is -1.35. Comparison of the MF of the IRAS 19343+2 026 region with tha t of the 



Trapezium cluster (slope = -1.21 for M/M0 > 0.6) measured by lMuench et al.l (120021 ) reveals 
a resemblance, in the sense that both MFs rise monotonically beyond M/M0 > 0.6. The 
cluster region slope is shallower than the Salpeter value (-1.35). This indicates that the 
star formation process is not yet complete in the region. Star formation is ongoing, although 
previous episodes of star formation activity have resulted in a significantly low mass content. 
This situation represents active star formation over a period of 3 - 5 Myr as observed in 
other examples of massive star forming regions (lAscenso et al.l 120071 ). 

We have also estimated the mass of the cluster by integrating the observed mass function 
above the completeness limit {logMn^ ~ 0.1 M0 from Fig. 7). A mass of ~ 307 Mq is 
obtained. For th e masses below the completeness limit, we have extrapolated the theoretical 
IMF slopes from iKroupal (J200ll ) down to the brown dwarf limit (see Fig. 7). The integration 
over the IMF yields a total cluster mass, Mtotai ~ 585 Mq. This is com parable to the masse s 



of clusters such as those associated with the Trapezium and Mon R2 (jLada &: Ladall2003l ) 



3.3 Spectroscopy 

Long-slit HK band spectra of two bright and five faint sources were obtained to investigate 
the nature of the stellar sources in the cluster. In Fig. 1, the straight lines mark the posi- 
tioning of the slit which, for each exposure, covered multiple stars. The numbers by the side 
of the line identify the stars seen through the slit. In Fig. 8 the spectra of each of these stars 
are shown. The stars 1, 2 and 3 are the bright stars that are close to the FIR emission peak. 
Stars 4, 5, 6 and 7 are fainter (presumably low mass) stars that lie in the central regions of 
the cluster. Stars 3, 5, 6 and 7 occupy the T Tauri zone in our JHK CC diagram in Fig. 4, 
while star 1 lies in the HAeBe zone of this plot. 

The two bright stars (1 & 3) that coincide with the FIR peak (see Fig. 1) display the HI 
Bracket series of lines indicating ionised emission. Note that star 2, which lies in between 
stars 1 and 3 on the slit, does not show the same recombination lines. The HI lines are 
thought to arise very close to the star, since the 120" long slit encompasses the dense nebula 
and we see no signature of extended nebular emission. As will be shown in the next section, 
the SED modelling of these bright objects classifies them as massive young stars. Star 1, the 
best massive protostellar candidate in the region, shows intense Paa emission, in addition 
to the HI recombination lines. Both stars, 1 and 3, display a rising slope long-ward of 2 
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Figure 7. The mass function (with itviV error bars) of the statisti cally cleaned IRAS 19343+ 2026 cluster region derived from 
the UKIDSS data. The plot is based on the evolutionary models of lLeieune fc Schaererl 1120011 ) and |Palla fc Stahler (1999 ), for 
an assumed age of 3 Myr. The vertical dashed a nd dot-dashed l ines represent 90% and 100% completeness limits, respectively. 
The dotted line shows the theoretical IMF from iKroupal 1120011) . for the masses below the completeness limit (logMi^-^ ~ 0.1 

M©). 



/im. Together, these features indicate that stars 1 and 3 are early B type stars capable of 
producing an ionised sphere around themselves (jOsterbrockl Il989l ) . Main-sequence B stars 
are characterised by HI absorption features which become weaker for more massive stars. 



However, for O type stars, he 



increasing mass (iHanson et al. 



ium a bsorption features appear and become prominent with 



I2OO5I). Neither HI nor He absorption features are visible in the 
bright stars 1 and 3, suggesting that the photospheres are likely obscured by the surrounding 
compact HII region and dense material. 
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In contrast, the remaining spectra of fainter stars (4, 5, 6 & 7) are featureless, and display 
a slope that is either flat or rising short-ward of 2 /im. The absence of any emission features 
suggests a relatively evolved T Tauri star (Antoniuc ci et al. 2008), • while the absence of 



any absorption features would suggest heavy veiling (iGreene &: Ladalll996l and references 
therein). Therefore, the fainter stars 4, 5, 6 & 7 are best matched by T Tauri type stars. The 
stars 6 & 7 are separated by a projected angle of 1.5"; both show a common spectral feature 
at 2.4 /im. This broad feature was not identified with any known spectral lines. The feature 
is likely to be a poorly removed telluric absorption feature (owing to poor sky beyond 2.3 
/im) and/or a local artifact in the slit. We note that there is no such spectral feature known 
in previous literature, even in the most embedded sources with lots of emission lines. 



3.4 SED modelling of bright stars 

The four bright stars that dominate the central region of the infrared nebula, namely stars 
1, 3, 8 and 9 (see Fig. 1), are resolved in the Spitzer IRAC bands for which point source 
photometry is available. Note that these four stars are the brightest NIR sources associated 
with the IRAS 19343+2026 cluster (see Table 1). Only stars 1, 3 and 8 were detected in the 
MIPS 24 /im images. Stars 3, 8 and 9 are also detected in our optical BVRI images. Of 
these, star 8 is a multiple star with unreliable photometry since the PSF is larger than the 
separation of the multiple components. Therefore, the optical magnitudes of star 8 are used 
as upper limits in the modelling. 

The SEDs of the s e four sources were modelled using the online SED fitting tool developed 



by iRobitaille et al.l ( 120071 ). The results of the SED modelling are summarised in Table 2. 
In Fig. 9, the observed data points (shown as black dots and filled triangles) and the fitted 
SED models are plotted. The solid black line shows the best fit model while the grey lines 
represent models that satisfy the criteria of x^ ~ X^best < 3, where x^ is per data point. The 
weighted means and standard deviations in Table 2 are calculated using all of the parameters 
from the models that satisfy the above criteria. The weight is the inverse of the x^ value. The 
standard deviations quoted in Table 2 help the reader to visualize the spread in the range 
of values between multiple models that satisfy the above criteria. The dotted line represents 
the photospheric emission assumed in the model. 

It is interesting to note that the FIR peak (star 1) is a young massive star of ~ 7.6 
Mq which is deeply embedded in an envelope. The remaining bright stars are 6-7 Mq objects 
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Figure 8. HK band spectra of the stars marked in Fig. 1. The shaded area represents the region of poor atmospheric 
transmission. 
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Table 1. Photometry used for SED modelling" 



RA (2000) 


DEC (2000) 


B 


V 


R 


/ 


J 


H 


K 


3.6/im 


4.5/^m 


5.8fi 


8.0^m 


deg 


deg 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


294.12900 


20.55662 


- 


- 


- 


- 


14.99 


13.18 


11.11 


8.10 


6.90 


5.68 


4.27 


294.12899 


20.55319 


- 


19.60 


17.50 


16.20 


12.68 


11.07 


9.98 


8.37 


7.70 


6.96 


- 


294.12634 


20.55140 


16.50 


15.00 


13.88 


13.22 


11.22 


10.60 


10.11 


9.25 


8.81 


8.31 


7.51 


294.12489 


20.55042 


20.50 


18.30 


16.70 


15.60 


12.09 


10.85 


9.86 


8.54 


8.11 


7.82 


7.36 



"0.1 mag (10%) error is assumed on all magnitudes. See text for details. 
''The UFTI J, H, K photometry of all the cluster members is available in the online table. 



with nearly revealed photospheres surrounded by remnant dust envelopes (indicated by the 
second peak in the SED models). These results are consistent with the fact that star 1 is not 
optically visible, unlike stars 3, 8 and 9. The VLA NVSS survey shows centimeter free-free 
continuum emission coincident with the central region of this nebula. The cm continuuum 
emission may represent the ionised gas in the central nebula associated with these four early 
B type stars. 



Star 1 


10.0 


7.6±0.6 


5.2±0.6 


-1.5±0.7 


1.2±0.7 


Star 3 


28.8 


6.3±0.8 


5.4±0.3 


-2.0±1.0 


0.8±0.7 


Star 8 


14.6 


6.1±0.6 


4.6±0.3 


-1.8±0.7 


0.9±0.6 


Star 9 


16.7 


7.4±1.8 


5.5±0.4 


-2.5±2.2 


-1.0±2.5 
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Table 2. SED modelling results 

Source x^ Mass (Mq) Log(Age) yr log(Mdisk) Mq log(Me„„) Mg log{Maisk) Mgyr"! log{Meuv) Mgyr-l 

-6.8±1.1 -4.3±0.6 

-7.2±1.2 -5.1±0.9 

-7.0±0.9 -3.8±0.4 

-7.2±2.3 -4.8±0.5 

4 CONCLUDING REMARKS 

In the previous sections we characterised the stellar content associated with the massive 
protostellar object candidate, IRAS 19343+2026. The stellar content is composed of at least 
four young (~ 10^ yr) B type stars and a rich population of low mass stars at ~ 1- 3 Myr. The 
FIR peak is modelled as a 7.6 Mq massive young star. The cluster is surrounded by a bright 
infrared nebula seen from 3-24 /im, indicative of dust and polycyclic aromatic hydrocarbons 
(PAH) emission. While the Spitzer IRAC and MIPS bands display the infrared nebula that 
becomes increasingly brighter from 3.6 yum to 24 /im, the NIR bands predominantly display 
the dense stellar population embedded in this nebula. The following facts can explain why 
the cluster population is better detected in deep K band images than in Spitzer images. 
The contribution of stellar photospheres rapidly falls off longward of 2 /im. Further, the 
background emission due to PAHs is strong in the bands longward of 4.5 /im. The above two 
facts, together with the relatively large PSF of the Spitzer data, reduces the sensitivity and 
contrast of the revealed low mass stars. However, the massive stars and the most embedded 
objects which appear bright at longer wavelengths are best revealed by the Spitzer images. 
This can explain why deeper NIR i^-band data such as the one presented here can unveil 
the lower mass population better than the GLIMPSE survey. 

The fainter and relatively older low mass population revealed by the NIR images are 
uniformly arranged over the infrared nebula (see Fig. 1), coinciding well with the warm 
dust distribution seen in 24 /im emission. As shown in Section 3.1, this population is best 
described to be 1 - 3 Myr old. In contrast, the four bright sources are found to be in the age 
range of 10^ - 10^ yr (see Table 2). This suggests that low mass star formation occured in the 



cluster prior to the formation of massiv e stars. Although this idea 



from observations of other regions (e.g. 



las be en argued before, 



Kumar. Tafalla &: Bachiller 



20041 ). we arrive at the 



same conclusion based on a more rigorous analysis using a CM diagram, HK-hand spectra, 
and SED modeling. If all embedded clusters are born with a universal IMF, then sampling 
the IMF at any given time, uniformly over all mass ranges, is naturally expected to result in 
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fewer massive stars and a greater number of low mass stars. Consequently, in the evolution 
of an embedded cluster (expected to follow an universal IMF), there exists an early phase 
when the total number of massive stars will be zero with a non-zero population of low mass 
stars. This simple statistical reasoning suggests that massive stars should appear after the 
low mass stars in a young cluster. Thus, it appears quite clear that the sequence for star 
formation is "low mass stars first and massive stars next", at least for clusters forming B 
type stars. The situation for clusters forming O type stars requires further investigation. 
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deg deg mag mag mag 
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15.11 


294.11145 


20.56095 
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16.05 
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16.56 
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17.20 


16.81 
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APPENDIX A: UFTI JHK PHOTOMETRIC DATA 

Only a few lines of Table Al are printed here. The complete form of the table is available 
online only. 



